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Abstract 

We compute the vector meson decay constants /d* from the simulation of twisted 
mass QCD on the lattice with Af = 2 dynamical quarks. When combining these 
values with the pseudoscalar £)( s )-meson decay constants, we were able (i) to show 
that the heavy quark spin symmetry breaking effects with the charm quark are large, 
Id* /fD s — 1-26(3), and (ii) to check the factorization approximation in a few specific 
5-meson non-leptonic decay modes. Besides our main results, fr>* = 278 ± 13 ± 
10 MeV, and fo* = 311 ± 9 MeV, other phenomenologically interesting results of 
this paper are: f D *JfD* = 1-16 ± 0.02 ± 0.06, /d;//d = 1-46 ± 0.05 ± 0.06, and 
f D Jf D * = 0.89±0.02±0.03. Finally, we correct the value for B(B° -> D+tt-) quoted 
by PDG, and find B(B° -> D+-R-) = (7.8 ± 1.4) x 10" 7 . Alternatively, by using the 
ratios discussed in this paper, we obtain B(B° ->■ D + ir~) = (8.3 ± 1.0 ± 0.8) x 10 - '. 
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1 Introduction 



Vector meson decay constants are important ingredients in the particle physics phenomeno- 
logical description of various processes pp. They are particularly handy when checking on 
the validity of the factorization approximation in non-leptonic decay channels involving a 
vector meson. Of particular interest are the 5°-meson non-leptonic modes in which one 
of the two charged mesons in the final state is D^l and the other one is a light meson. 
Factorization is expected to work quite well for this class of decay modes [2]. In the spe- 
cific cases such as B° — y Di*} tt + , the factorization was shown to be exact in the limit of 
infinitely heavy quark mass [3] . Away from that limit and in other similar situations, such 
as B° — y D^ + 7T~, the factorization is an assumption [6]. While it is difficult to check on 
the extent to which the factorization approximation works for the absolute values of the 
branching fractions, it is relatively simple to do it when considering the ratios of various 
modes. i?-factory experiments at BaBar and Belle provided us with many accurate mea- 
surements of the non-leptonic 5-decays, of which particularly interesting are the measured 
B° — y Z}|*j + 7r~, and B° — y D^) + D~ modes jl]. For example, in the naive factorization the 
amplitude for the B° —y D + n~ decay writes, 



where the Wilson coefficients ci^/i) are known to next-to-leading order (NLO) in QCD 
perturbation theory. The hadronic matrix elements, instead, are non-perturbative quanti- 
ties and need to be accurately evaluated on the lattice. The amplitudes for the other modes 
can be written similarly, and it is easy to see that in the suitable ratios the semileptonic 
hadronic matrix element |(7r~|67^u|.B )| cancels out, leaving only the ratios of vacuum-to- 
D^fl meson matrix elements to be computed j6]. This is a favorable situation for the lattice 
QCD studies that will be explored in this paper. The factorization approximation is shown 
to work for B° — y Z}j*j~ 7r + . We will present the results of our computation of the D* and 
D* meson decay constants (fo* and /dj), that can be combined with the pseudoscalar 
ones (fo and fo s ) and check the extent to which the factorization approximation provides 
the adequate description of the ratios of several non-leptonic B decay modes. From the 
lattice QCD point of view, particularly appealing is the computation of Jd*/ Id s because 
both the physical charm and strange quarks are directly accessible in the simulations, and 
the final physical result requires only a smooth extrapolation in the sea light quark mass 
and to the continuum limit. 

Our computation of fr>* and fn* is made on the ensembles of gauge field configurations 
produced by the European Twisted Mass Collaboration (ETMC) at four different lattice 
spacings and for several light sea quark masses [TJ, by using the twisted mass QCD on the 
lattice (LtmQCD) [8]. In this way the chiral and continuum extrapolations are well under 
control. The results presented here are unquenched, with iVf = 2 flavors of dynamical mass- 
degenerate light quarks. We briefly discuss the corresponding pseudoscalar meson decay 
constants as well that will help us answering to the question about the size of the heavy 

^ore specifically, from ref. [5], one reads: cf LO (ra 6 ) = -0.285(14) and c^ LO (m b ) = 1.132(8). 
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quark spin symmetry breaking effect in the heavy-light mesons when the heavy quark is 
charm. In our work this means that we want to measure how far from the static limit we 
actually are. In the static limit, lim (/d*//d) = 1. As we shall see, our results suggest 

that this ratio is considerably larger than one. 

The remainder of this paper is organized as follows. In Sec. |2] we remind the reader 
of the definition of the pseudoscalar and vector meson decay constants, and discuss the 
peculiarity of their lattice determination in the framework of LtmQCD. In Sec. [3] we sum- 
marize the details of the lattices used in our analysis, present our results as extracted from 
the correlation functions studied on each of the available data-sets, and after the chiral 
and continuum extrapolation we present our physical results. Using our results, we then 
compare the ratios of several non-leptonic decay modes measured experimentally with the 
results of the factorization approximation. In Sec. E]we briefly conclude. 

2 Open charm mesons and their decay constants 

The decay constants of charmed pseudo-scalar and vector mesons are defined through the 
following matrix elements, 

(0|c(0)7 M 7 5 g(0)|£>,(p)) = / D ,p MJ 

(2) 

(0|c(0) 7 ^(0)|D*(p,A)> = f D *m D *e^ 

where the index q stands for either the strange or the light w/d-quark. We assume isospin 
symmetry and do not distinguish between the u- and <i-quark mass. In the above expres- 
sions, p and are the meson momentum and the vector meson polarization, respectively. 

Instead of the axial current matrix element, in LtmQCD, it is far more convenient to 
extract the pseudoscalar meson decay constant from the matrix element of the pseudoscalar 
density. At the maximal twist, the renormalization constant of the quark mass and the 
pseudoscalar densities cancel exactly, and therefore 

(fi q + fic)(0\c(0h^mD q ( P )) = f Dq m 2 Dq , (3) 

with fi q and \x c being the light and charm bare quark masses respectively. The results of the 
computations of fu in LtmQCD at the maximal twist were first reported in ref. [9]. They 
were later corroborated in ref. [TO], i.e. after including the results of simulations made at 
finer lattice spacing, corresponding to /3 = 4.2. In this paper we will focus onto the vector 
meson decay constants, fo* , the results of which we present together with the pseudoscalar 
ones, as to better emphasize the fact that the charm is indeed far from the static heavy 
quark limit. In contrast to the pseudoscalar case, the vector meson decay constant cannot 
be extracted without explicitly accounting for a renormalization factor [Hj. The relevant 
constant for the physical vector current in LtmQCD is Z and its value for all our 
lattices has been computed non-perturbatively in ref. [TT] . 

The access to the vector meson mass (m D *) and its decay constant (/d*) on the lattice 
is made through the study of the large time separation in the two-point correlation function 
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involving the vector current, V$ = Z^ig^ciiQ, namely 



Cvvit) = (£Vi(£; t)V?(0) 0)> = -Tr S c (0, 0; x, t)jiS q (x, t; 0, 0) Ti 



X 



t > cosh[m J; (T/2-t)] 
m D * 



(0|V i (0)|D*(0,A)> e 



2 



m D *T/2 



(4) 



In the above expression S q ^(x; 0) stands for the light (charm) quark propagator. Recently, 
the authors of ref. [13] noted that the use of stochastic source propagators in the computa- 
tion of vector meson properties is not as advantageous as it is with the pseudoscalar mesons. 
That was in fact already observed in refs. [T^l[T5j . but in ref. [13] the authors show that the 
extraction of the light vector meson properties from the lattice can be improved if several 
stochastic sources in the propagator inversion are used. 

3 Lattice details and results 

This section contains the main details of our computation, including the discussion of the 
extrapolation procedure that leads to the phenomenologically relevant results. 

3.1 Lattices and correlation functions 

Details about the sets of gauge field configurations with N{ — 2 dynamical quarks that are 
used in this work are summarized in table [TJ 

Like in the case of the light-light vector currents, in the corresponding heavy-light 
correlation functions we observe that the point source propagators lead to the equally good 
results for the decay constants and hadron masses as those obtained by using the stochastic 
source propagators. We did not attempt using more stochastic sources. We merely note 
that the quality of the signal for the two-point correlation functions obtained by two kinds 
of propagators is the same and the effective mass plateaus remain indistinguishable, c.f. 
fig. [TJ □ The effective mass is the solution of, 



In the following the results we present are obtained by using the stochastic propagators [7]. 

Another worry when dealing with heavy quarks on the lattice is to make sure that the 
lowest lying state has been well isolated at moderately large time separations between the 

2 The behavior of my (t) obtained with the point source propagators is less stable even if the error bars 
are smaller, while the time dependence of the effective mass is more stable from the correlation functions 
with stochastic propagators but with slightly larger errors. As a result the vector meson masses from the 
correlators with either kind of propagators are indistinguishable. 
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3.8 


3.9 


3.9 


4.05 


4.2 


4.2 


L 3 x T 


24 3 x 48 


24 3 x 48 


32 3 x 64 


32 3 x 64 


32 3 x 64 


48 3 x 96 


# meas. 


240 


240* 


240 


240 


240 


96 


Mseal 


0.0080 


0.0040 


0.0030 


0.0030 


0.0065 


0.0020 


/4ea2 


0.0110 


0.0064 


0.0040 


0.0060 






Msea3 




0.0085 




0.0080 






Msea4 




0.0100 










a [fm] 


0.098(3) 


0.085(3) 


0.085(3) 


0.067(2) 


0.054(1) 


0.054(1) 


^a(5o 2 ) PHE21 


0.746(11) 


0.746(6) 


0.746(6) 


0.772(6) 


0.780(6) 


0.780(6) 


Ms 


0.0194(12) 


0.0177(11) 


0.0177(11) 


0.0154(10) 


0.0129(10) 


0.0129(10) 


Mc 


0.2331(82) 


0.2150(75) 


0.2150(75) 


0.1849(65) 


0.1566(55) 


0.1566(55) 



Table 1: Summary of the details about the lattice ensembles used in this work. The asterisk in the number 
of gauge field configurations used from the simulations at j3 = 3.9 is there to indicate that the data-set 
with /iscai = 0.0040 contains 480 configurations, while the sets with other values of /z se a contain 240. Data 
obtained at different /3's are rescaled by using the Sommer parameter ro/a, and the overall lattice spacing 
is fixed by matching /„■ obtained on the lattice with its physical value, leading to ro = 0.440(12) fm (c.f. 
ref. [12]). Strange and charm bare quark masses, fj, a and fi c respectively, are obtained as discussed in 
ref. [12]. All quark masses are given in lattice units. 

interpolating operators. In refs. [HI [10] the values for rrio q and fr> were obtained from the 
correlation functions with the local sources. We checked that by using Gaussian smearing 
(see e.g. ref. [IB]) the plateaus indeed become larger but the error bars on the resulting 
decay constants remain the same. Furthermore the central values for the pseudoscalar 
masses and decay constants remain unchanged (within the statistical accuracy). A similar 
situation holds true in the case of vector mesons, i.e. the effective mass indeed exhibits 
plateau 'earlier' when the smeared source operators are used. At moderately large time 
separations, where the extraction of the masses and decay constants is made, the effective 
masses obtained from correlation functions with local and with smeared interpolating oper- 
ators coincide. This point is illustrated in fig. [2] In what follows we will quote as our main 
results those obtained by using local operators with stochastic sources. We checked that 
our final results remain unaltered when the smeared correlation functions are combined 
with the local ones. 
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Figure 1 : Comparison of the effective masses of the heavy-light vector mesons obtained from the correla- 
tion function (Tj| computed by combining either point source propagators or stochastic source propagators. 
Illustration is provided for the data sets obtained at f3 = 3.9 and with fi sca = 0.0064. 



3.2 f D * and f D */fD s 

With the above comments in mind, we now present our results for the charm-strange pseu- 
doscalar and vector mesons. Since the ensembles of gauge field configurations used in this 
work are obtained with Nf — 2 light quarks, the results presented in this subsection are 
only partially unquenched. The strange valence or charmed quarks are directly accessible 
from our lattices but they do not have their "sea" -quark counterparts. Only the light sea 
quarks (denoted in what follows as m q ) have been included in the QCD vacuum fluctua- 
tions. The argument that the strange quark mass is too heavy to significantly contribute 
to the QCD vacuum fluctuations seems to be confirmed by the fact that the fu B decay 
constants obtained from the simulations with Nf = 2 (THl E] and Nf = 2 + 1 [18] are 
essentially equal. Of course such a comparison is obscured by the fact that different lattice 
QCD actions have been used in obtaining the respective results. As far as the charm quark 
is concerned, none of the currently available lattice results on charm physics include the 
charm quark in the sea. That is likely to change in the near future as the ETMC and MILC 
are continuously producing the gauge field configurations with Nf = 2 + 1 + 1 dynamical 
quark flavors, i.e. with both the strange and charm quarks in the sea [T^l I2U| . 

In tab. [2] we give the values for the masses and decay constants as obtained in each 
of the ensembles of gauge field configurations described in tab. HJ We should stress that 
the charm and strange quark mass were obtained from completely separated studies [T2] . 
and that the values of vector meson masses and decay constants are the net prediction 
of lattice QCD. For the reader's convenience we converted our results from lattice units 
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Figure 2: Effective mass of the heavy- light vector meson, extracted from the two-point correlation 
function ([4J obtained with either the local vector currents or by implementing Gaussian smearing. At 
larger time separations we see that the two results are totally compatible. Illustration is provided on the 
data at j3 = 3.9 and /i sca = 0.0100. 



to the physical ones by using the lattice spacing values given in tab. HJ From tab. [2] we 
also see that the errors on the vector meson decay constants are twice larger than for the 
pseudoscalar ones. Finally, to get to the physical result we need to extrapolate in the sea 
quark from the masses used in the simulations down to the physical w/rf-quark mass in 
the continuum limit. To do so we combine all results obtained with the valence physical 
charm and strange quarks, and fit our data to 

$(L > :) latt. = (/,, v /77777) !> '^- [1 + A s m^ + B s a 2 } , (6) 

where on the left hand side we combine the masses and decay constants into the quantity 
$(-D*) = f'D*y/^D*- A s and B s are the fit parameters. One can use the same formula 
to extrapolate the decay constant alone. An expression similar to eq. (|6]) is also used for 
the pseudoscalar decay constant, and so one can estimate the size of the heavy quark spin 
symmetry breaking effects in the decay constants through the ratio <&(£)*) /$(.D S ). Clearly, 
from our data, shown in fig. [3j we do not see any sea quark mass dependence. Furthermore 
the 0(a 2 ) effects are very small, and we finally obtain 

f D * = 311 ± 9 MeV , ^ = 1.26 ±0.03. (7) 

JD a 

Our result in eq. (J7J) shows that the heavy quark spin symmetry breaking effects are 
larger than 20% for the case in which the heavy quark is charm. This simply means 
that for the decay constants the charm physics is far away from the static limit in which 
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L 


/^sea 




m D . 


fD s 




3.80 


21 


/^seal 


1.773(54) 


2.008(59) 


0.280(8) 


0.352(15) 






/^sea2 


1.771(54) 


2.005(58) 


0.279(9) 


0.351(15) 


3.90 


24 


/^seal 


1.806(44) 


2.023(46) 


0.265(7) 


0.326(10) 






A*sea2 


1.807(44) 


2.029(49) 


0.268(7) 


0.336(15) 






/^sca3 


1.798(44) 


2.020(49) 


0.263(6) 


0.330(10) 






/^sea4 


1.810(44) 


2.029(49) 


0.272(7) 


0.333(15) 


3.90 


32 




1.803(44) 


2.022(48) 


0.266(6) 


0.331(10) 






/4ea2 


1.803(44) 


2.019(49) 


0.267(6) 


0.330(10) 


4.05 


32 


/^seal 


1.877(27) 


2.099(31) 


0.264(5) 


0.341(9) 






/^sea2 


1.878(27) 


2.079(32) 


0.267(5) 


0.331(9) 






/^sea3 


1.880(27) 


2.102(30) 


0.267(6) 


0.339(9) 


4.20 


32 


/^seal 


1.909(38) 


2.104(39) 


0.266(5) 


0.320(10) 


4.20 


18 


/^seal 


1.888(39) 


2.096(38) 


0.251(5) 


0.316(8) 



Table 2: Charmed meson decay constants and hadron masses with the light quark fixed to the strange 
quark mass [12]. All results are given in GcV. 

lim (/d*//d s ) — 1) an d that power corrections, proportional to 1/m™, are important. To 

m c — ¥oo s 

better appreciate this fact we remind the reader that f p + / f w + ~ 1.58, and f K *+ / f K + ~ 1.36 
(see for example (UETJ). Note also that both fo* and fn* / /d s are somewhat larger than the 
results obtained in the quenched approximation [22]. We should also add that from a similar 
chiral and continuum extrapolation of the meson mass we obtain mc. = 2141(22) MeV, 
compatible with the experimentally established m^,? G = 2112 MeV [I]. 

3.3 f D *, Id* /Id, /d;//d*, and other ratios 

Next we present our results for masses and decay constants of charmed mesons containing 
a light (up or down) quark. We treat the light quark mass as fully unquenched, i.e. we keep 
the light valence and sea quarks degenerate in mass. As before, and for the readers' conve- 
nience, we list all our results in physical units and together with the corresponding masses 
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Figure 3: Dependence of the ratio of $>(D*)/$(D S ) on the sea quark mass, where 
The lattice discretization errors appear to be negligible. 



= mf s (2 GeV). 



and decay constants of the pseudoscalar mesons. As in the previous subsection, we proceed 
by combining all our data to make the combined chiral and continuum extrapolation via 
the expression, 



latt. 



\Phys 



1-^ 



4(4vr/ 7r ) 2 



m: 



\ogml + A u4 ml + B Utd a 



where again we considered the quantity /d* ^ m D* = ®(D*) for which the above formula, 
derived in heavy meson chiral perturbation theory (HMChPT), can be used to guide the 
chiral extrapolation [23]. Since this step turns out to be the dominant source of system- 
atic uncertainties, we should spend a few more lines to discuss it. Despite the fact that 
the window of light quark masses used in our data stretches down to about 1/7 of the 
strange quark mass, the fact that the chiral logarithmic correction to this quantity is large 
introduces an important uncertainty in the result of the extrapolation. First of all, when 
applying the above expression to the data obtained with the charmed heavy quark, it is 
not clear what value for the coupling of the heavy-light mesons to a soft pion, g, to use. 
While its value in the static heavy quark limit (m c — > oo) has been computed on the lattice 
and was shown to be g ~ 0.4 0.5 [21], its value computed with the propagating charm 
quark was found to be g ~ 0.7 [25] . Since the chiral corrections were computed in the 
static heavy quark limit one is tempted to use the value for g obtained in the static limit 
as well. However, knowing that the l/m™-corrections are sizable, the validity of HMChPT 
in the charm sector is dubious and one should try other possibilities, including the larger 
g ~ 0.7, or even the simple linear chiral extrapolation. Another aspect, that should not 



L 






m D * 


fn q 


fo* 


3.80 24 


/^seal 


1.722(52) 


1.946(58) 


0.264(9) 


0.321(19) 




/^sea2 


1.735(52) 


1.967(57) 


0.267(9) 


0.338(19) 


3.90 24 


/^seal 


1.741(44) 


1.969(48) 


0.244(8) 


0.311(21) 




A*sea2 


1.751(43) 


1.970(51) 


0.251(7) 


0.313(21) 




/^sea3 


1.748(43) 


1.967(49) 


0.248(6) 


0.311(11) 




/^sea4 


1.774(42) 


1.992(51) 


0.262(7) 


0.317(21) 


3.90 32 


/^seal 


1.733(43) 


1.958(51) 


0.244(7) 


0.312(20) 




/4ea2 


1.732(42) 


1.935(48) 


0.244(7) 


0.298(11) 


4.05 32 


/^seal 


1.808(27) 


2.024(31) 


0.241(7) 


0.312(12) 




/4ea2 


1.823(27) 


2.026(33) 


0.250(6) 


0.313(11) 




/^sea3 


1.840(27) 


2.067(30) 


0.256(8) 


0.328(10) 


4.20 32 


/^seal 


1.866(35) 


2.064(38) 


0.252(6) 


0.306(12) 


4.20 48 


/^seal 


1.806(37) 


2.032(42) 


0.224(7) 


0.290(15) 



Table 3: Charmed non-strange meson masses and decay constants for various light quark masses. All 
results arc obtained with the light valence quarks mass degenerate with the sea quark mass. Like in tab. [21 
all numbers are given in physical units [GcV] . 

be ignored, is the fact that the HMChPT formulas were obtained by taking into account 
the lowest doublet of heavy-light mesons only. As we know, however, the first orbital ex- 
citations (scalar and axial mesons) turned out to be very close to the lowest lying states 
(pseudoscalar and vector mesons), and their inclusion might modify the chiral extrapola- 
tions as well [26J . To make a fair assessment of the chiral uncertainties, each of our central 
values will correspond to the result obtained by using the expression (jSJ) with g = 0.45, 
and the difference between that and: (a) the result obtained by using g = 0.7 in eq. (jSJ), 
(b) the result obtained through the linear chiral extrapolation, will be added separately as 
an uncertainty to our result. To give the reader a good feeling about the smoothness of 
the chiral and continuum extrapolations, here are the numbers obtained from the fit to our 

$(£>*) : fn*V^ = 382(2) MeV 3 /\A u4 = 0.0(3) GeV" 2 , B u4 = 0.0(3) fm" 2 , (9) 
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which resemble very much the shape of the combined chiral and continuum extrapolations 
in the case of pseudoscalar mesons, namely, 



HD q ) : f Dv ^B = 276(8) MeV 3 / 2 , A' u4 = 0.0(2) GeV" 2 , B' u4 = 0.0(2) fm" 2 . (10) 

It goes without saying that the chiral logarithmic correction derived in HMChPT in the 
static limit is the same for both the pseudoscalar and vector mesons. 

What might look striking from the above numbers is the fact that the discretization 
errors are essentially absent. This is actually a fortunate artifact of the quantities $(£**) 
and $(-D g ). The results listed in tab. [3] obviously show the non-negligible a 2 -effects. The 
decay constant exhibits a positive slope in a 2 and the continuum value for f D M is smaller 
than those obtained at fixed value of the lattice spacing. More pronounced slope in a 2 
shows up in the results for the heavy-light meson masses, but of the sign opposite to that 
in the decay constants. Those two effects cancel against each other to a large extent, 
which is why = f D (*) y/m D (*) is manifestly insensitive to discretization effects (the 

parameter B^ d obtained from our fit is consistent with zero.) 
Our final results are: 

fo* = 275 ± 13^7 MeV , ^ = 1.28 ±0.06. (11) 

3d 

where the second error corresponds to the uncertainty arising from the chiral extrapolation: 
the upper error is the difference with the central value of the linear extrapolation, and the 
lower is the difference with result obtained by using g = 0.7 instead of g = 0.45 in eq. ([8]). 
As for the D* meson mass, we obtain mp* = 2041(22) MeV, that is again somewhat larger 
than the measured m™ G = 2010 MeV. We should mention that the previous studies with 
tmQCD on the lattice observed that the light vector meson masses are larger than the 
physical ones [T4"| IT5]. Although to a much lesser extent, this feature seems to persist with 
the heavy-light mesons too. [f| 

Other interesting and phenomenologically useful ratios are: 

p- = 1.17 ±0.02+°^, ^ = 1.19 ±0.01+^, 

3D* 3D 

fp- = 0.90±0.02+° 2 , fp- = 1.48±0.05±°|. (12) 

3D* 3D 

Note, in particular, that our results are consistent with the recently claimed inequality 
Id* I Id* < IdJ/d deduced from HMChPT, after including the power corrections and 
assuming the validity of SU(3) chiral perturbation theory for the light quarks 



3 For example in ref. [14] the authors quote mx* = 0.981(33) GeV, the central value of which is about 
10% larger than m^- 1 ? 6 = 0.892 GeV. In our case, the central value for our mo* is 2% larger then the 
experimentally measured mass. This result is obtained by using the smeared source operators. The non- 
smeared (local) sources give a value compatible with the smeared one, within the error bars. Its central 
value is however a tad larger and is 3.5% higher than the experimentally established £>*-meson mass. 
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Figure 4: SU(3) breaking ratio for the charmed vector meson decay constants. Dashed curves correspond 
to the chiral extrapolation using cq. ((SJ) for each lattice spacing and their spread indicates the size of 
discretization effects. The full thick curve is the result of chiral extrapolation in the continuum limit. 



4 Verifying the factorization approximation 

With the information on the charmed meson decay constants in hands, we can now illustrate 
on the specific examples how one can check the validity of the factorization approximation 
in several particular non-leptonic 5-decays mentioned in introduction of the present paper. 
To begin with, we remind the reader that the following modes have been measured [4]: 



B(B° D+Ti-] 
B(B° -»• D* + n-] 

B(B° -> D+7T-] 
B(B° D* s + n-] 



(4.6 ±0.4) x lCr 5 , 
N.A., 

(2.16 ±0.26) x 10~ 5 , 
(2.1 ±0.4) x 10~ 5 , 



B(B° -»■ D + D~) = (2.11 ± 0.31) x 10~ 4 
B(B° -> D* + D~) = (6.1 ± 1.5) x 10" 4 
B(B° ->■ D+D-) = (7.2 ± 0.8) x 10~ 3 
B(B° -»■ D* s + D-) = (7.4 ± 1.6) x 10" 3 



(13) 



We consider the following four ratios: 



Ri 



Bo. 



B(B° £>+ 



7T 



5(5° D+7T 

5(5° -> D+D- 
B(B° -> D+D- 



Vcs_ 

V cd 

Vrs 



V, 



cd 



\(m B ,m Da ,m n ) 
_ \{m B ,m D ,m w ) 
\(m B ,m Ds ,m D ) 
\(m B ,m D ,m D ) 



1/2 



1/2 



^0 



(m 2 



rn 



D a 



fDs 

fD 

fDs 

fD 
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R i 



B(B° ->■ Z>+tt 

B(B° -> D*+D- 
B(B° -> D*+D- 



cd 



\(m B ,m D *,m n 

\(m B ,m D *,m D ) 
\(m B ,m D *,m D )_ 



3/2 



m 



D* 



1 2 



rn 



m 



D* 



fD* 



(14) 



where on the right hand side we explicitly write the expressions obtained in the factoriza- 
tion approximation, with the usual X(a, b, c) = [a 2 — {b — c) 2 ][a 2 — {b + c) 2 ]. Ratios of the 
semileptonic form factors on the right hand side can be taken to be one to excellent accu- 
racy, so that by using our results from eq. (JTj), together with V cs /V c d = 1/ tan 6^ = 4.32, 
we obtain: 



^(fact.) 
p (fact.) 



R, 
R 



(fact.) 
3 

(fact.) 



26.0 ±0.4 ±2.6, vs. 

25.7 ±0.4 ±2.6, vs. 

23.8 ±0.8 ±2.5, vs. 
22.7 ±0.8 ±2.4, vs. 



^(exp.) 
r»(exp.) 



R 
R 



(exp.) 
3 

(exp.) 



0.47 ±0.07. 

34.1 ±6.3, 
N.A. , 
12.1 ±4.0. 



(15) 



The first ratio shows a very large disagreement between the factorization approximation 
and experiment, and it is mainly due to B° — > D + ir~ for which the experimental branch- 
ing fraction is by almost two orders of magnitude larger than the one obtained in the 
factorization approximation. To understand the origin of this discrepancy we looked 
in the original paper by the Belle Collaboration [28], and realized that the value for 
B(B° — > D + tt~) quoted by PDG, and given in eq. ( |T3|) above, is not actually measured 
but deduced from the measurement of B(B° — > D+ir~), and by imposing the validity of 
factorization, similar to our ratio Ri above. Furthermore, PDG erroneously identified the 
measured value of R Dw = 1.71(11)(9)(2)% [23], with R Dt = B(B° D + n-)/B(B° 
D~7i + ), while in fact the definition given in ref. [28], after assuming the factorization, 
reads R DlT = ^/B(B° D+n-)/B(B° D--K+). We corrected that error, and by using 
B(B° D-n + ) = (2.68 ± 0.13) x 10" 3 [1], we obtain 



B(B° D- 



7T 



\PDG 



(7.8 ± 1.4) x 10' 



(16) 



Since the value of B(B° — > D + n ) is not measured independently but extracted by impos- 
ing the factorization approximation, the comparison of our ^ fact -) with corrected i?[ exp ' ) 
would not provide us with any useful information. Instead, we can follow the same recipe 
as PDG, and use our value for ^ fact -) from eq. ( !T5|) . to arrive at 



B(B° D~ 



7T 



L3±1.0±0.8) x 10- 



(17) 



where the first error is the experimental uncertainty coming from B(B° .D+tt"), and the 
second is our theoretical uncertainty within the factorization approximation. The reason 
why the results in eqs. (fT6|) and (jT7|) differ is the fact that in ref. [28] the value of SU(3) 
breaking ratio f B J f B = 1.164(11) has been used, while our value is foj Id = 1.18(1) (6). 
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Furthermore, in converting our flf** to the result (02) we used B(B° — > D+ir ) from 
ref. jl], which is the average of the results obtained by Belle [2S] and BaBar [29J . 

Similarly, our result for the ratio Rf 1101 ^ can be combined with the measured B(B° — > 
D* + ix~) to deduce, 



B(B -»■ D* + ir' 



(8.8 ± 1.6 ±0.9) x 10' 



where the first error reflects the experimental uncertainty in B(B° — > D* + ir~), and the 
second is theoretical error within the factorization approximation. 

In the case with two charmed mesons in the final state, we see that the factorization 



approximation works rather well. The difference between and R^^'' can be used 

to constrain the corrections to the factorization coming from either the penguin (Cabibbo 
suppressed) contributions or the final state interaction [30] . 

To use other results from eqs. (17ffT2|) and check whether or not the factorization approx- 
imation in ratios of various decay modes agrees with experiment, one also needs different 
form factors, and an extra assumption is needed. For example, 



,(exp.) 



RtL 



B(B° D*+tc- 
B(B° D+7T 



1.09 ±0.05) 



fact. 



[\{m B ,m D * 



m 7 



3/2 



(m 2 B 



m 



2\2 



[A(mB,m D ,,m^)] 



1/2 



[m 



m 



D s 



m 



Dt 



7 B^ { 



m 



f. 



D 3 



(19) 



where, in the last line, we used our result from eq. ([7]). Knowing that the two semileptonic 
B — > ir form factors are slowly varying functions at low q 2, s, and that F + (0) = Fq(0), it is 
reasonable to assume that F^ n (m 2 r 



F B_Mr ir 



2 D ), and therefore: 



R 



(fact.) 



1.09 ±0.05, vs. R 



(exp.) 



1.0 ±0.2 



(20) 



a comparison that again goes in favor of the validity of factorization. 



5 Summary 

In this paper we discussed the extraction of the charmed vector meson decay constants from 
the simulations of tmQCD on the lattice with iVf = 2 dynamical light quarks. The results 
we obtain are presented along with those of the pseudoscalar mesons, clearly showing that 
the heavy quark spin symmetry breaking effects are large in the case of the heavy charm 
quark, i.e. that corrections oc 1/m" are sizable. The dominant sources of errors are those 
coming from the chiral extrapolations when the extrapolation in the valence light quark 
mass was needed (i.e. for D-mesons). We do not make a 'guestimate' of the systematic 
uncertainty that might arise from the non-included strange and charm quark in the sea. 
That point will be numerically assessed once the ongoing simulations with A^f = 2 + 1 + 1 
dynamical quark flavors are completed. As for the main numerical results presented in this 
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paper, after symmetrizing the error bars, we quote: 



fry, = 278 ± 13 ± 10 MeV , f 



311 ± 9 MeV 



1.26 ±0.03 




1.28 ±0.06, 



1.16 ±0.02 ±0.06. (21) 



In some cases the results obtained in this paper can be used to verify the validity of the 
factorization approximation in non-leptonic B-decay modes. We also corrected the error in 
PDG and instead of the currently reported B(B° — > D + ir~) = (4.6 ± 0.4) x 10~ 5 , we obtain 
B(B° — > D + 7i~) = (7.8 ± 1.4) x 10~ 7 , which is in very good agreement with the value we 
obtain from another ratio [Ri in eq. ffT4l)]. namely B(B° — > D + n~) = (8.3 ± 1.0 ± 0.8) x 
10~ 7 . Both values are obtained by relying on the validity of the factorization approximation. 
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